9 Fluid dynamics and Rayleigh-Benard convection 



In these lectures we derive (mostly) the equations of viscous fluid dynamics. 
We then show how they may be generalized to the problem of Rayleigh- 
Benard convection — the problem of a fluid heated from below. Later we show 
how the RB problem itself may be reduced to the famous Lorenz equations. 

The highlights of these lectures are as follows: 

• Navier-Stokes equations of fluid dynamics (mass and momentum conser- 
vation) . 

• Reynolds number 

• Phenomenology of RB convection 

• Rayleigh number 

• Equations of RB convection 

Thus far we have dealt almost exclusively with the temporal behavior of a 
few variables. 

In these lectures we digress, and discuss the evolution of a continuum. 

9.1 The concept of a continuum 

Real fluids are made of atoms or molecules. 

The mean free path fmfp is the characteristic length scale between molecular 
collisions. 

Let I/hydro be the characteristic length scale of macroscopic motions. 
Fluids may be regarded as continuous fields if 

hydro -^mpf- 
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When this condition holds, the evolution of the macroscopic field may be 
described by continuum mechanics, i.e., partial differential equations. 



To make this idea clearer, consider a thought experiment in which we measure 
the density of a fluid over a length scale £ using some particularly sensitive 
device. We then move the device in the x-direction over a distance of roughly 

m. 



Suppose i ^ Li ^ £mpf- Then we expect the density to vary greatly in space 
as in Figure (a) below: 



e 




(b) 




hydro 



We expect that the fluctuations in (a) should decrease as £ increases, (statistics 

tells us that these fluctuations should decrease like where N (x is the average number of 

molecules in a box of size i. ) 

On the other hand, if ^ ~ ^hydro (see (c)), variations in density should reflect 
density changes due to macroscopic motions (e.g., a rising hot plume), not 
merely statistical fluctuations. 

Our assumption of a continuum implies that there is an intermediate scale, 
£ ~ L2, over which fluctuations are small. Thus the continuum hypothesis 
implies a separation of scales between the molecular scale, Li ~ f^fp, and the 
hydrodynamic scale, Lhydro- 

Thus, rather than dealing with the motion ~10^^ molecules and therefore 
~6 X 10^^ ordinary differential equations of motion (3 equations each for 
position and momentum), we model the fluid as a continuum. 



83 



The motion of the continuum is expressed by partial differential equations for 
evolution of conserved quantities. We begin with the conservation of mass. 



9.2 Mass conservation 



Let 



p = density 
u = velocity 



of a macroscopic fluid particle 



Consider a volume V of fluid, fixed in space: 




ds is an element of the surface, |ds| is its area, and it points in the outward 
normal direction. 

u is the velocity. 

The outward mass flux through the element ds is 

pu ■ ds. 

Therefore, 

rate of mass loss from V = pu ■ ds. 

J s 

The total mass in V is 

/ pdv 
Jv 

Thus the rate of mass loss may be rewritten as 



Shrinking the volume, we eliminate the integrals and obtain 

J pu ■ ds/V 



dp y 

— = - hm 

at v^-^o 
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Recall that the RHS above is the definition of the divergence operator. We 
thus obtain 

We see that to conserve mass, a net divergence creates a corresponding change 
in density. 

For incompressible fluids, 

p ~ constant. 

(This result is not an assumption, but instead derives from the assumption that the Mach number, 
the square of the ratio of the fluid velocity to the speed of sound, is much less than unity.) 

Then 

V • w = 0. 

which is the equation of continuity for incompressible fluids. 



9.3 Momentum conservation 

We seek an expression of Newton's second law: 

-^(momentum of fluid particle) = force acting on fluid particle (22) 

yJiL/ 

9.3.1 Substantial derivative 



We first focus on the LHS of (22). 

There is a conceptual problem: ^(particle momentum) cannot be given at a 
fixed location, because 



• the momentum field itself changes with respect to time; and 

• fiuid particle can change its momentum by flowing to a place where the 
velocity is different. 
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To better understand this problem physically, consider how a scalar property — 
the temperature T — of a fluid particle changes in time. 

A small change 5T is produced by a small changes 5t in time and 5x^ 6y, 6z 
in the position of the fluid particle: 

dT^^ dT^ dT^ dT^ 
6T = —6t + —6x + —6y + —6z 
ot ox oy oz 

Divide by 5t to obtain the rate of change: 

5T dT OTSx OTSy dTSz 
— = \ \ ^ -\ 

6t dt dx 6t dy 6t dz 6t 

In the limit 6t ^ 0, 

6x 6y Sz 

¥ ^ ^ ^ 

The rate of change of T of a fluid particle is then 



DT dT dT dT dT 

+ Uj:— + Uy— + U, 

+ U-VT 



Dt dt ' dx ' ^ dy ^ dz 
dT 



dt 
where 

D^ = a^ + ^-^ 

is the substantial derivative or convective derivative operator. 

Thus we see that the temperature of a fluid particle can change because 

• the temperature field changes "in place" (via d/dt); and 

• the particle can flow to a position where the temperature is different (via 

V). 

Note that the same analysis applies to vector fields such as the velocity u: 

Du du , ^ ^ 
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Therefore the velocity u enters Dw/Dt in 2 ways: 



• u changes (in place) as the fluid moves (d/dt) 

• u governs how fast that change occurs {u • V) . 

This dual role of velocity is the essential nonlinearity of fluid dynamics and 
thus the cause of turbulent instabilities. 

We can now express the rate-of-change of momentum per unit volume (i.e., 
LHS of (22)): 

Du du ,^ ^ 

p is outside the differential because a fluid particle does not lose mass. Density changes thus mean 
volume changes, which arc irrelevant to the momentum change of that particle. Above we have 
written the (rate of change of momentum) per unit volume, which need not be equal to the rate of 
change of (momentum per unit volume). 



9.3.2 Forces on fluid particle 



To obtain the full dynamical equation, we need the RHS of 

p—— = Force acting on fluid particle / unit volume. 



These forces are 

• body force (i.e., gravity) 

• pressure 

• viscous friction (internal stresses) 

— * 

Body force. We represent the externally imposed body force by F. 
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Pressure. Fluid flows from high to low pressure. Thus 

pressure force dp 
unit volume dx 



in 1-D 



= -Vp 



in 3-D 



Viscous friction. Viscous stresses are the source of dissipation in fluids. 
They resist relative movements between fluid particles. 

For example, the shear flow 



is resisted more by high viscosity fluids than low viscosity fluids. 

This resistance derives from molecular motions. (A nice analog is Rcif's picture of 

two mail trains, one initially fast and the other initially slow, that trade mailbags.) 

In the simple shear flow above, there is a flux of a:-momentum in the y- 
direction. 



In Newtonian fluids^ this flux, which we call Pxy^ is proportional to the gra- 
dient: 

_ _ dux_ 

where 77 is called the dynamic viscosity, r] has units of mass/ (length x time). 

The shear stress can occur at any orientation. Analogous to the 1-D Newto- 
nian condition above, we deflne the viscous momentum flux 



Pij = -V 



duj 
dxj 
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The conservation of momentum requires that the divergence of the momen- 
tum flux Pij be balanced by a change in the momentum of a fluid particle. 
Loosely stated, 



2 



dt 

We thus find that 



VISCOUS j ^-^3 j ^-^j 



viscous force 

= rjV u. 



unit volume 

(a careful derivation requires consideration of the tensorial relationship between viscous stress and 
the rate of deformation.) 

Newton's second law then gives the Navier-Stokes equation for incompressible 
fluids: 

^ ^ , 

(mass per unit vol) x acceleration 



stresses on fluid element per unit vol body force per unit vol 



Incompressibility arose from our negelect of compressive forces on fluid ele- 
ments. 



9.4 Nondimensionalization of Navier-Stokes equations 

Deflne the characteristic length scale L and velocity scale U. We obtain the 
non-dimensional quantities 

I X t y I ^ 
"=f7' *=*!• 

The dynamical equations (without body force) become 

du' _ ^ 1 ^/2 -^ 



^^, + (fl'.v'K = -vV + -v 
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where 

Re = Reynolds number = 
is the dimensionless control parameter. 

The Reynolds number quantifies the relative importance of the nonlinear term 
to the viscous term. To see why, note the following dimensional quantities: 

\pu ■ Vu| ~ — — nonlinearity 
hV^wl ~ dissipation 



Their ratio is 



\pu-Vu\ pUL o 11 . 

= Reynolds number 



High Re is associated with turbulence (i.e., nonlinearities). Low Re is asso- 
ciated with laminar or creeping flows dominated by viscous friction. 

Note that as long as Re remains the same, the dimensional parameters like 
U and L can change but the the flow (i.e., the equation it solves) does not. 
This is dynamical similarity. 



An example is running vs. swimming: 

= 0.01 cm^/sec 



= 0.15 cm^/sec and ' 



air 



water 



On the other hand, comparing 100 meter world records, 

10"^ cm o , 

f^run ~ -r^ = 10 cm/sec 

10 sec 

10^ cm ^ / 

L/swim ~ -73 ~ 2 X 10 cm/sec 

55 sec 

Taking L ~ 100 cm, 

Re(swim) ~ 2 x 10^ and Re(run) ~ 6 x 10' 
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Thus for both swimming and running, Re ~ 10^, well into the turbulent 
regime. Surprisingly, despite the slower speed of swimming, Re(swim) is 
somewhat greater. 

Another example: bacteria swimming in water is roughly like us swimming in 
molasses, since the the small size and slow speed of bacteria would correspond 
to a larger and faster body in a more viscous fluid. 



9.5 Rayleigh-Benard convection 

In a thermally expansive fluid, hot fluid rises. 

R-B convection concerns the study of the instabilities caused by rising hot 
fluid and falling cold fluid. 

Typically,, fluid is confined between two horizontal, heat-conducting plates: 



8 



T=Tq (cold) 



A 

d fluid 

\l_ 



T=Tq+ hT (hot) 




pure 

conduction 



temperature 



In the absence of convection — the transport of hot fluid up and cold fluid 
down — the temperature gradient is constant. 

Two cases of interest: 

• 5T small: no convective motion, due to stabilizing effects of viscous 
friction. 



• 5T large: convective motion occurs. 
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How large is a "large ST" ? We seek a non-dimensional formulation. 
The following fluid properties are important: 

• viscosity 

• density 

• thermal expansivity 

• thermal diffusivity (heat conductivity) 

Convection is also determined by 

• (i, the box size 

• 6T (of course) 

Consider a small displacement of a cold blob downwards and a hot blob 
upwards: 



T=Tg+ 6r 

Left undisturbed, buoyancy forces would allow the hot blob to continue rising 
and cold blob to continue falling. 

There are however damping (dissipation) mechanisms: 

• diffusion of heat 

• viscous friction 
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Let Dt = thermal diffusivity, which has units 

length 



[Dt] = 



time 



The temperature difference between the two blobs can therefore be main- 
tained at a characteristic time scale 



We also seek a characteristic time scale for buoyant displacement over the 
length scale d. 



Let 



po = niean density 

Ap = — apoAT, a = expansion coefficient 
Setting AT = 6T, 

buoyancy force density = \gAp\ 

= gapo6T. 

Note units: 

r rn-n maSS 

[gapodT\ = 



(length) 2 (time) 2 

The buoyancy force is resisted by viscous friction between the two blobs 
separated by ~ d. 

The viscous friction between the two blobs diminishes like 1/d (since viscous 
stresses oc velocity gradients) . The rescaled viscosity has units 



mass 



(length) 2 (time) 



Dividing the rescaled viscosity by the buoyancy force, we obtain the charac- 
teristic time Tm for convective motion: 

-q/d Tj 



buoyancy force go-PQ d 6T 
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Convection (sustained motion) occurs if 

time for motion < diffusion time for temperature difference 

Thus convection requires 

— > constant 

or 

— — — 5T = Ra > constant 
'qUT 

Ra is the Rayleigh number. A detailed stabiUty calculation reveals that the 
critical constant is 1708. 

Our derivation of the Rayleigh number shows that the convective instability 
is favored by 

• large 6T, a, d, pq. 

• small 77, Dt- 

In other words, convection occurs when the buoyancy force p^gaS 6T exceeds 
the dissipative effects of viscous drag and heat diffusion. 

Note that box height enters Ra as d^. This means that small increases in box 
size can have a dramatic effect on Ra. 



9.6 Rayleigh-Benard equations 
9.6.1 Dimensional form 

We employ the Boussinesq approximation: density perturbations affect only 
the gravitational force. 

The momentum equation is therefore the Navier-Stokes equation augmented 
by the buoyancy force: 

^ + w • W = -— Vp + uV'^u - ga(T - To) 
ot po 
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Here we have written the kinematic viscosity 

z/ = ri/po 

The mass conservation equation is again 

V-u = 0. 

We now additionally require an equation for the convection and diffusion of 
heat: 

dT 

— + {u-V)T = DtV^T. 

9.6.2 Dimensionless equations 

The equations are nondimensionalized using 

length scale = d 
time scale = (f /Dt 
temperature scale = (5T/Ra. 

An additional dimensionless parameter arises: 

Pr = Prandtl number = z^/-Dt, 
which is like the ratio of momentum diffusion to thermal diffusion. 
We shall employ the dimensionless temperature fluctuation 

6 = deviation of dimensionless T from the simple conductive gradient 

The mass conservation equation is 

V-w = 

Momentum conservation yields (i is a unit upward normal) 

= - Vp + + V^w 



1 

Pr 



du -L^ 
at 
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The heat equation becomes 

39 

— j^U-Ve = Ra(u • z) + V^^ 
Note that there are two nonhnear terms: 

— » 

• u • Vu 

• u-ve 

Their relative importance depends on Pr: 

— » 

• small Ft ^ u • Vu dominates. Instabilities are "hydrodynamic." 

• large Ft ^ u • V9 dominates. Instabilities are thermally induced. 

9.6.3 Bifurcation diagram 

For Ra < Rac, there is no convection. 

For Ra > Rac, but not too large, a regular structure of convection "rolls" 
forms, with hot fluid rising and cold fluid falling: 

T=To 

r=ro+6 —d — 

Now imagine placing a probe that measures the vertical component v of 
velocity, somewhere in the box midway between the top and bottom. A plot 
of t'(Ra) looks like 
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V 




+ 



convection (stable) 







conduction 



Ra 



rest 



V 



Such a plot is called a bifurcation diagram. Here the stable states are bold 
and the unstable states are dashed. 

Note that we cannot know in advance whether the velocity will be up or 
down. This is called symmetry breaking. 

9.6.4 Pattern formation 

Rayleigh-Benard convection makes fascinating patterns. Some examples: 

• Figures 22.3-8, Tritton. 

• Plate 1, Schuster, showing quasiperiodic regime. (The 40 sec period is 
not precise: note details in upper right are not quite periodic.) 

• Plumes: Figure 22.12, Tritton 

• Plumes in the wind: Zocchi, Moses, and Libchaber (1990) 

• Collective plumes: Zhang et al (1997). 

9.6.5 Convection in the Earth 

The Earth's radius is about 6378 km. It is layered, with the main divisions 
being the inner core, outer core, mantle, and crust. 
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The Earth's crust — the outermost layer — is about 30 km thick. 
The mantle ranges from about 30-2900 km. 

The mantle is widely thought to be in a state of thermal convection. The 
source of heat is thought to be the radioactive decay of isotopes of uranium, 
thorium, and potassium. Another heat source is related to the heat deriv- 
ing from the gravitational energy dissipated by the formation of the Earth 
roughly 4.5 Ga. 

At long time scales mantle rock is thought to flow like a fluid. However its 
effective viscosity is the subject of much debate. 

One might naively think that the huge viscosity would make the Rayleigh 
number quite small. Recall, however, that Ra scales like d^, where d is the 
"box size". For the mantle, d is nearly 3000 km!!! 

Consequently Ra is probably quite high. Current estimates suggest that 

3 X 10^ < Ra,,antle < lO' 

which corresponds to roughly 

10^ X Rac < Ra^antie < lO^Rac 

The uncertainty derives principally from the viscosity, and its presumed vari- 
ation by a factor of about 300 with depth. 

Some pictures illustrate these ideas: 

• Science cover, 26 May 1989 

• van der Hilst seismic tomography, showing cold slab descending toward 
the core-mantle boundary. 

• Gurnis (1988) simulation/cartoon showing breakup of continents. 

• Zhang and Libchaber (2000) showing floater-plates. 



Thermal convection is the "engine" that drives plate tectonics and volcanism. 
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It turns out that volcanism is, over the long-term, responsible for the CO2 
in the atmosphere, and thus the source of carbon that is fixed by plants. 
(Weathering reactions remove C from the atmosphere.) 

Thus in some sense thermal convection may be said to also sustain life. 

That is, without convection, there probably would be no CO2 in the atmo- 
sphere, and therefore we wouldn't be around to discuss it... 
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